Current constraints on the epoch of cosmic acceleration 
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The cosmographic expansion history of the universe is investigated by using the 557 type la 
supernovae from the Union2 Compilation set along with the current estimates involving the product 
of the CMB acoustic scale £a and the BAO peak at two different redshifts. Using a well-behaved 
parameterization for the deceleration parameter, q(z) — qo+qiz/(l + z), we estimate the accelerating 
redshift z acc = —qo/(qo (at which the universe switches from deceleration to acceleration) and 

investigate the influence of a non-vanishing spatial curvature on these estimates. We also use the 
asymptotic value of q(z) at high-z to place more restrictive bounds on the model parameters qo and 
qi, which results in a more precise determination of the epoch of cosmic acceleration. 

PACS numbers: 98.80.-k, 95.36.+X, 98.80.Es 
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Introduction- Although a late-time cosmic acceleration 
is becoming observationally well established (at > 12a 
from some independent analyses [l| ) , the physical mech- 
anism behind this phenomenon constitutes a completely 
open question nowadays in Cosmology. In principle, this 
accelerating expansion may be the result of unknown 
physical processes involving either the existence of new 
fields in high energy physics, the so-called dark energy 
in the context of general relativistic models @, or the 
infrared modification of gravity one should expect, e.g., 
from extra dimensional physics, which would lead to a 
modification of the effective Friedmann equation at late 
times 3, 4] . Another possibility of modification of gravity 
is to add terms proportional to powers of the Ricci scalar 
R to the Einstein-Hilbert Lagrangian, the so-called f(R) 
gravity models (see, e.g., Q and references therein). 

Given the current state of uncertainty that remains 
over which mechanism may provide a complete descrip- 
tion for cosmic acceleration, a better understanding of 
some characteristics of this phenomenon can certainly 
help us distinguish among many alternative models of 
universe. This is the case, for instance, of the duration 
of the accelerating phase, which depends crucially on the 
physical mechanism of acceleration. In principle, it ma" 
distinguish between thawing and freezing potentials 
or between some classes of braneworld models and the 
standard ACDM scenario Q , and may also test the com- 
patibility between classes of dark energy models and the- 
oretical constraints from String/M theory 1 . 

Another important aspect in this discussion concerns 
the epoch of the transition from an initially decelerated 
to an accelerating phase. A purely kinematic record of 
the expansion history of the universe, without regard to 



1 As discussed in Ref. [7| , an eternally accelerating universe seems 
not to be in agreement with String/M-theory predictions, since 
it is endowed with a cosmological event horizon which prevents 
the construction of a conventional S-matrix describing particle 
interactions. This is the case of freezing scenarios as well as of 
our current standard ACDM model. 



its cause or the validity of any particular metric theory 
of gravity, can provide answers to basic questions con- 
cerning the acceleration history and model-independent 
constraints with which one may test cosmological mod- 
els. As an example, we note that one of the features of 
evolving dark energy is that the redshift at which the uni- 
verse switches from deceleration to acceleration, z acc , is 
remarkably different from that in ACDM scenario with 
the same amount of dark energy today, a fact that is 
manifested in the Cosmic Microwave Background (CMB) 
data as a modified integrated Sachs- Wolfe effect [8(. 

Our goal in this Brief Report is to derive updated cos- 
mographic bounds on the epoch of the cosmic accelera- 
tion, without making any assumption on the energy con- 
tent of the Universe (for similar analyses using either 
different parameterizations and expansion formulas for 
the deceleration parameter q(z) or different data sets, 
see [l|, H, 0)- Our approach consists on a parametric 
approximation of the deceleration parameter along the 
cosmic evolution, given by q(z) = q + q\z/{\ + z) [l|. 
In order to ensure a period of structure formation dur- 
ing the matter-dominated era, we also make use of the 
asymptotic value of q(z) at high redshift to constrain the 
model parameters qo and q\. The statistical analysis is 
performed using the most recent SNe la observational 
data, the so-called Union2 sample of 557 events [Tjj, to- 
gether with the recent estimates involving the product 
of the CMB acoustic scale I a and the baryonic acoustic 
oscillation (BAO) peak [ll|. In what follows, we outline 
the main assumptions for our analysis and discuss our 
main results. 

Basic equations and Analysis-We start our discussion 
by assuming the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) line element. In such a background, the comov- 
ing distance to a given redshift z can be written as 
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where X(u) = exp {-/ "[l + g(it)]<21n(l + u)} and the 
function Sk is denned by one of the following forms: 
Sk( r ) = sinh(r), r, and sin(r) for open, flat and closed 
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geometries, respectively. In the above equation, we 
have used the definition of the deceleration parameter 
q = —da /a 2 or, equivalent ly, 

, . 1 d\nH 2 

where ait) = (1 + z)^ 1 stands for the cosmological scale 
factor, H = a/a is the Hubble parameter and a dot de- 
notes derivatives with respect to time. 

In order to proceed further, we follow [l| and adopt in 
our analysis a two-parameter expression for the deceler- 
ation parameter 

q(z) = q + qi— — , (3) 
1 + z 

where go is the current value of the deceleration parame- 
ter and in the infinite past q(oo) = qo + q\. Clearly, this 
parametric form for q(z) was inspired from one of the 
most popular parametrization of the dark energy equa- 
tion of state [12| and, although very simple, seems to be 
flexible enough to mimic the q(z) behavior of wide class of 
accelerating models [PH. FromEq. ©, we can also cal- 
culate the Hubble parameter for the above parametriza- 
tion, i.e., 



H(z) = H (l + z) 1+qo+q 




The accelerating redshift (q(z acc ) — 0) for parameteriza- 
tion ^ is given by 



^acc — , ■ V*-v 

qo + q\ 

In order to derive constraints on the epoch of cosmic 
acceleration, we use in our analyses the most recent SNe 
la compilation available, namely, the Union2 sample of 
Ref. yi|- This compilation is an update of the origi- 
nal Union sample that comprises 557 data points in the 
redshift range 0.015 < z < 1.4 including recent large 
samples from other surveys and using SALT2 for SN la 
light-curve fitting. 

The predicted distance modulus for a supernova at red- 
shift z, given a set of parameters s, is /j, p (z\s) = m — M = 
51ogdz, + 25, where m and M are, respectively, the ap- 
parent and absolute magnitudes and cLl stands for the 
luminosity distance, dz, = (1 + z)r(z) (in units of mega- 
parsecs). In our analysis, we estimated the best fit to the 
set of parameters s = (ffoj<7o,9i) by using a % 2 statis- 
tics Xsnc = Yh=i [^ P i z \ s ) -/40)] /°f) where ^(z\s) 
is given above, (J? a (z) is the extinction corrected distance 
modulus for a given SNe la at Zi, and Oi is the uncer- 
tainty in the individual distance moduli. The Hubble 
parameter H D is considered a nuisance parameter so that 
we m arg inalize over it by using the analytical method of 
Ref. [ijj. 

Along with the SNe la data, and to impose more re- 
strictive bounds on the parameters qo and q%, we use 



the two major inputs involving acoustic oscillations that 
come from the CMB and baryon oscillations data. We 
follow Ref. [ll| and use constraints derived from the 
product of the CMB acoustic scale I a = ndA(z*)/r s (z*) 
and the measurement of the ratio of the sound hori- 
zon scale at the drag epoch to the BAO dilation scale, 
T s (zd)/ Dv(zbao)- In the above expressions, d,A(z«) is 
the comoving angular-diameter distance to recombina- 
tion z* = 1089 and r s (z*) is the comoving sound hori- 
zon at photon decoupling given by r s (z*) = j^jdz, 
which depends upon the speed of sound before recom- 
bination, c s . Zd — 1020 is the redshift of the drag 
epoch (at which the acoustic oscillations are frozen 
in) and the so-called dilation scale, Dy, is given by 
Dy{z) — [czr 2 (z)/H(z)] 1 ^ 3 . By combining the ratio 
r s (z d = 1020)/r s (z* = 1090) = 1.044 ±0.019 Q with 
the measurements of r s (zd) / Dv(zbao) at zbao == 0.20 
and 0.35 from Ref. [l5j], Sollerman et al. fTT| found 

/0.20 = d A (z*)/D v (0.2) = 17.55 ± 0.65 

/o.35 = d A (z*)/D v (0.35) = 10.10 ± 0.38 . 

In order to derive the constraints on acceleration red- 
shift z acc discussed in the next section we minimize the 
function Xt = XsNc + Xcmb/bao> where Xcmb/bao = 

[/o.2(z|s) - /o. 2 ] 2 /<Jq.2 + Lf0.35O|s) - /o.35] 2 /CTo.35- 

Results-The results of our analysis are shown in Figs. 
1 and 2 and in Table 1. Figure la shows 68.3% and 
95.4% confidence regions in the (go, qi) plane for param- 
eterization §3§ arising from Union2 SNe la sample (blue 
ellipses) and measurements of CMB/BAO ratio (green 
bands) when vanishing spatial curvature (flka — 0) is as- 
sumed. Note that, although not very restrictive when 
applied separately, these SNe la and CMB/BAO con- 
tours constrain different regions in the parametric space, 
showing the complementarity of these data sets, and re- 
sulting in very tight bounds on the (50, 3i) space when 
used together. Quantitatively, if we use only SNe la data 
we obtain z acc = 0.75±0.35 (2a) whereas the CMB/BAO 
measurements provide z acc — 1.20 ±0.10 (2a). From the 
SNe la + CMB/BAO analysis, we find z acc = 0.87 ±0.10 
(2a) with reduced xt — 0.97, which is roughly consistent 
with the ACDM prediction with fi m = 0.27 ± 0.04 [3, 
i.e., z acc € [0.64,0.88]. The constraining power of the 
joint SNe la + CMB/BAO analysis is clearly shown by 
the 2a black contours. In Fig. lb, we show the evolution 
of the deceleration parameter q(z) for the resulting 2a 
intervals of go and gi obtained using the data sets con- 
sidered above. We note that the general behavior of q(z) 
is very similar to the one predicted by several classes of 
accelerating scenarios and reflects the flexibility of pa- 
rameterization ([3]) discussed earlier. For comparison, the 
ACDM curve with Sl m = 0.27 is also shown (thick line). 

For most of the viable cosmological scenarios, we ex- 
pect the Universe to be matter-dominated at early times 
(i.e., after the radiation dominance), which implies q — 
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FIG. 1: a) Confidence regions in the (go, qi) plane for pa- 
rameterization ([3j) arising from Union2 SNe la sample (blue 
ellipses) and measurements of CMB/BAO ratio (green bands) 
for flko = 0. Black regions correspond to the resulting bounds 
from the joint analysis. Contours are drawn for A\ 2 = 2.30 
and Ax 2 = 6.17. b) Evolution of the deceleration parameter 
q(z) for the resulting 2a intervals of qo and q\ obtained using 
the data sets discussed in the text. 

1/2. Thus, in order to ensure this past dark matter- 
dominated epoch (and the standard a(t) oc t 2 ^ 3 law), 
whose existence is fundamental for the structure for- 
mation process to take place, we assume the constraint 
q(z » 1) = 1/2. In this case, q + qi = 1/2 and Eq. ^ 
becomes 

?(*) = + I) ^ , (6) 

which has the advantage of reducing our analysis to one- 
parameter fitting. From this equation, the acceleration 
redshift is now simply given by z acc = — 2qo. Note also 
that, for z » 1, Eq. (gj) reduces to H oc (1 + z) 3 / 2 , as 
expected for a decelerating matter-dominated universe. 

Figure 2a shows the behavior of the deceleration pa- 
rameter as a function of z when the above constraint 
is taken into account. Note that, differently from the 
results shown in Panel lb, the CMB/BAO estimate on 
z acc is now moved to a much lower value, i.e., z acc = 
0.32 ± 0.20 (2a) whereas the Union2 sample alone gives 
Zacc — 1.16 ± 0.18 (2a). The joint analysis involving 
these two data sets plus the above constraint on q(z) at 
high-z limit provides z acc — 0.71 ± 0.12 (2a). This is 
in good agreement with the concordance ACDM model 
prediction with VL m = 0.27 (z acc = 0.75) as well as 
with a wCDM model 2 with the same amount of non- 
relativistic matter today and equation-of-state (EoS) pa- 
rameter ranging in the interval —1.76 < w < —0.63. 



2 The accelerating redshift in tcCDM models with a constant EoS 
parameter w can be obtained from the expression: Sl m (l + 
Zacc) 3 + (1 + 3w)(l -fl TO )(l + z acc j i{1+w) = 0, which reduces to 
the ACDM prediction z acc = (2r2 A /Sl m ) 1 / 3 - 1 when to = -1. 
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FIG. 2: a) The same as in Fig. lb when the high-z limit 
on the deceleration parameter q(z >> 1) — 1/2 is consid- 
ered, b) Confidence regions in the (go, flko) space from SNe 
la + CMB/BAO + high- z limit on q(z). Note that a spa- 
tially closed and current accelerating universe is favored by 
the joint analysis (black contours). All the contours are drawn 
for Ax 2 = 2.30 and A X 2 = 6.17. 



Further, by considering the above high-z limit of q(z), 
we show in Fig. 2b the effect of a non-vanishing spatial 
curvature on the z acc estimates. From this analysis, we 
see that SNe la and CMB /B AO data provide almost or- 
thogonal bounds on the (go,^fco) plane and, therefore, 
very tight constraints on the current values of the de- 
celeration and curvature parameters. In particular, we 
found q = -0.61 ± 0.08 and Q k0 = -0.08 ± 0.02 at 2a 
level, which in turn favors a spatially closed and current 
accelerating universe, with z acc = 1.21±0.17 (2a). 

For the sake of completeness, we also derive the accel- 
erating redshift by assuming the current la bounds on 
the curvature parameter from WMAP7 results, flho = 

-0.00231^^6 03 • As can be seen from Table r > the 
effect of curvature for this realistic interval of ilko is 

small, being < 8% for the two extreme cases of open 
(flko = 0.0031) and closed (f2fco = —0.0079) geometries. 
The main quantitative results of this paper are summa- 
rized in Table 1. 

Conclusions-We have performed a cosmographic anal- 
ysis to estimate the epoch of cosmic acceleration with- 
out making any assumption about the underlying grav- 
itational theory and energy components of the universe. 
We have used in our analyses the most recent SNe la 
data, the so-called Union2 sample of 557 events in the 
range 0.015 < z < 1.4, along with estimates involving 
the product of the CMB acoustic scale I a and the BAO 
dilation scale (Dy) at two different redshifts (z = 0.2 and 
z = 0.35). 

Wc have shown that the complementarity of these two 
data sets (see Figs, la and 2b) poses tight constraints on 
Zacc (with uncertainties of ~ 15%), with the current esti- 
mates being in good agreement with the standard ACDM 
scenario with Q m ~ 0.27. In agreement with previous 
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TABLE I: Current limits on z acc from SNe la and CMB/BAO data. 



Test 


<7o 


gi 


Zacc 




Xmin 


n k0 = -0.0079 


CMB/BAO 


-1.36 ± 0.02(1ct) ± 0.04(2a) 


2.48 ± 0.03(la) ± 0.07(2a) 


1.21 ±0.05(1(7) 


± 0.11(2(7) 


1.18 x 10" 4 


Union2... 


. -0.66 ± 0.03(1ct) ± 0.07(2cr ) 


1.53 ±0.19(lcr) ±0.38(2cr) 


0.76 ±0.18(lcr) 


±0.36(2cr) 


543.22 


Union2 + CMB/BAO... 


. -0.63 ±0.01(lcr) ±0.03(2(7 ) 


1.34 ±0.02(lcr) ±0.05(2a) 


0.89 ± 0.04(lcr) 


±0.10(2cr) 


543.97 


Q k0 = 0.0031 


CMB/BAO 


-1.36 ± 0.02(1ct) ± 0.04(2a) 


2.51 ± 0.03(1(7) ±0.07(2a) 


1.18 ±0.05(la) 


±0.10(2(7) 


5.68 x 10" 5 


Union2 


-0.66 ±0.03(lcr) ±0.07(2a) 


1.55 ± 0.19(lcr) ±0.38(2a) 


0.74 ±0.17(lcr) 


± 0.34(2cr) 


543.22 


Union2 + CMB/BAO.... 


. -0.63 ±0.01(lcr) ±0.03(2cr) 


1.36 ± 0.02(lcr) ± 0.05(2a) 


0.86 ± 0.03(lcr) 


±0.10(2cr) 


543.98 


fi M =0 


CMB/BAO... 


. -1.39 ±0.02(lcr) ±0.04(2ct) 


2.55 ± 0.03(lcr) ± 0.07(2a) 


1.20 ± 0.05(la) 


±0.10(2cr) 1.51 x 10~ 3 


Union2... 


. -0.66 ± 0.03(lcr) ± 0.07(2a ) 


1.54 ± 0.19(lcr) ±0.38(2a) 


0.75 ±0.17(lcr) 


±0.35(2ct) 


543.22 


Union2 + CMB/BAO.... 


. -0.63 ±0.01(lcr) ±0.03(2cr) 


1.35 ± 0.02(lcr) ± 0.05(2a) 


0.87 ± 0.04(lcr ) 


±0.10(2cr) 


544.00 


O fc0 = and q(z » 1) = 1/2 


CMB/BAO 


-0.16 ±0.05(1ct) ±0.10(2cr) 




0.32 ±0.10(lcr) 


±0.20(2ct) 


1.76 


Union2 


-0.58 ±0.04(lcr) ±0.09(2cr) 




1.16 ±0.08(lcr) 


±0.18(2cr) 


543.79 


Union2 + CMB/BAO.... 


. -0.36 ±0.03(lcr) ±0.06(2a) 




0.71 ±0.06(lcr) 


±0.12(2cr) 


591.59 



analyses, our study strongly favors a Universe with a 
recent acceleration (q < 0) and an earlier decelerated 
phase (qi > 0). We have also checked the influence of 
curvature on the estimates of the deceleration parameter 
and accelerating redshift (see Table I). 

Finally, it is worth emphasizing that the analyses per- 
formed here clearly provide independent evidence for a 
dynamical model in which a deceleration/acceleration 
transition happened at z ~ 1. From the combination 
of data used in our analyses it is also possible to detect 



a cosmic acceleration today with great confidence, even 
allowing for arbitrary curvature. 
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